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SUMMARY 


The  concept  of  utilizing  optical  energy  (radiation  pressure)  to  provide 
motion  control  and  stability  for  water  droplets  in  the  size  range  applicable  to 
atmospheric  cloud  physics  (5-100  pm)  is  investigated.  It  is  shown  that  the  inher- 
ent stabilizing  forces  of  a gaussian  laser  beam  (TEMqo)  together  with  the  beam's 
ability  to  manipulate  individual  particles  makes  a powerful  tool  for  studying  drop- 
let or  ice  crystal  growth  and  for  initiating  relative  droplet  motion  in  order  to 
investigate  various  collision/coalescence  phenomena.  The  utility  of  optical  mo- 
tion control  (OMC)  will  be  greatly  enhanced  with  the  advent  of  Shuttle/Spacelab 
payloads  such  as  the  Atmospheric  Cloud  Physics  Laboratory  (ACPL)  in  which  the  am- 
bient opposing  accelerations  will  be  reduced  to  at  least  1/1000  of  that  encountered 
on  earth.  This  future  ACPL  capability  together  with  current  programs  for  develop- 
ing shuttle  compatible  lasers  in  the  100-200  mW  range  will  make  it  feasible  to  con- 
sider cloud  physics  experiments  employing  OMC  on  the  earliest  flights  of  the  ACPL. 

Specific  experimental  capability  of  OMC  in  the  Shuttle/Spacelab  environment 
will  include: 


\ 


1.  Providing  velocities  ranging  from  5 to  < 0.01  cm/sec. 

2.  Providing  precise  positioning  (±  5 pm)  of  water  droplets  up 
to  130  pm  radius  (droplet  growth  experiments) . 

3.  Providing  simulation  of  1 g conditions  (velocity  or  accel- 
eration) for  droplets  up  to  'v*  13  pm  radius  and  fractional 
g levels  for  larger  droplets. 

Representative  ACPL  experiments  which  can  utilize  these  capabilities  are  discussed. 

In  addition  to  the  ACPL  objective,  a program  for  developing  the  OMC  concept 
for  use  in  an  earth-based  laboratory  is  presented.  It  is  shown  that  OMC  is  appli- 
cable to  experiments  involving  the  coalescence  efficiencies  of  cloud  droplets  in 
the  growth  transition  region  between  7 and  15  pm  radius.  In  particular,  the  in- 
fluence of  droplet  charge,  electric  fields,  surfactants  and  turbulence  on  the 
collision/coalescence  process  in  this  size  range  can  be  studied.  The  proposed  ter- 
restrial laboratory  program  will  also  provide  a means  of  evaluating  and  implementing 
the  OMC  concept  for  experiments  in  the  Shuttle/Spacelab  ACPL  facility. 
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1.0 


BACKGROUND  AND  INTRODUCTION 


1.1 


Past  and  Current  Work 


The  growth  of  single  cloud  droplets  and  the  interaction  between  cloud 
droplets  has  long  been  recognized  as  an  area  of  investigation  which  must  be  pur- 
sued in  order  to  more  fully  understand  the  formation  of  precipitation  on  the 
microphysical  level  (Eaton,  1971).  The  fundamental  question  is  "How  do  the  small 
drops  (1-50  microns)  in  clouds  grow,  in  the  time  period  of  one  hour  or  less,  to 
the  size  of  rain  drops  (100-1000  microns)?  A NASA  sponsored  facility,  the  Atmo- 
spheric Cloud  Physics  Laboratory  (ACPI.) , is  currently  being  designed  with  the  ob- 
jective of  providing  answers  to  questions  of  this  nature. 

The  ACPL  is  a Shuttle/Spacelab  payload  consisting  of  a one  meter  wide 
double  rack  containing  three  chambers  along  with  the  equipment  required  to  pro- 
vide thermal  and  moisture  control,  as  well  as  aerosol  preparation,  conditioning 
and  characterization  capability.  Earth  orbit  will  provide  ambient  accelerations 
of  below  10~3  g (Final  Report  (ACPL),  1977). 


The  initial  objectives  are  directed  toward  relatively  simple  cloud  forma- 
tion experiments  involving  primarily  the  warm  cloud  nucleation  phenomena.  This 
experiment  selection  was  made  as  a result  of  taking  into  consideration  the  new- 
ness of  experimentation  in  a weightless  space  environment  which  in  turn  dicated 
the  selection  of  a relatively  simple  and  better  understood  experiment. 

The  ACPL  is  being  designed  to  permit  future  growth  and  expansion  of  the 
initial  experiment  cabability.  The  primary  areas  that  have  received  the  highest 
scientific  priority  are  those  dealing  with  scavenging  and  those  involving  ice. 
Some  dynamic  experiments  (e.g.,  collision/coalescence)  may  also  benefit  from  the 
low  acceleration  environment. 


Most  experiments  beyond  nucleation  will  require  the  observation  of  parti- 
cles greater  than  a few  micrometers  for  periods  of  time  greater  than  a hundred 
seconds.  The  dynamic  requirement  of  some  experiments  (i.e.,  relative  velocities 
between  droplets),  coupled  with  the  fact  that  the  ACPL  will  not  be  in  a true  zero 
acceleration  environment,  dictate  that  methods  will  be  required  for  position  and 
motion  control  of  droplets. 
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A number  of  motion  control  methods  are  available  and  many  of  them  are 
being  considered  for  "zero-gravity"  applications.  Each  technique  requires  spe- 
cific ambient  and  particle  characteristics,  for  example:  (a)  acoustical  methods 
require  a gas  (no-vacuum)  environment  and  particle  sizes  approaching  the  wave- 
length of  the  sound  (usually  above  one  centimeter),  (b)  electrostatics  and  elec- 
trodynamics require  charged  bodies,  (c)  electromagnetics  requires  specific  con- 
ductivity properties,  and  (d)  optical  methods  require  non-absorbing  particle 
properties.  The  first  three  of  the  above  operate  on  all  particles  within  a volume 
while  the  fourth,  optical,  can  be  preferentially  controlled  so  that  individual 
particles  can  be  uniquely  manipulated.  The  physical  properties  of  water  and  ice 
coupled  with  the  requirements  of  the  atmospheric  microphysics  experiments  (both 
particle  size  and  dynamics)  indicate  that  optical  motion  control  (OMC)  is  the  most 
suited  approach. 

The  concept  of  "radiation  pressure"  (not  to  be  confused  with  radiometric 
pressure)  was  first  introduced  by  Kepler  in  1619  when  he  suggested  that  the  pres- 
sure of  sunlight  caused  comet  tails  to  always  point  away  form  the  sun.  Until  the 
relatively  recent  invention  of  the  laser  in  1960,  radiation  pressure  was  not  con- 
sidered as  a force  to  be  reckoned  with  in  the  laboratory.  When  the  tremendous 
light  intensities  available  from  relatively  low  power  (^  0.1  W)  lasers  are  applied 
to  small  particles  ('v  10  pm  radius)  the  resulting  forces  are  competitive  with  the 
previously  dominant  ones  such  as  radiometric  forces  (due  to  thermal  gradients  in 
the  medium)  and  Earth's  gravity.  The  feasibility  of  utilizing  radiation  pressure 
for  controlling  the  movement  of  small  particles  is  further  enhanced  by  the  poten- 
tial for  conducting  the  experiments  in  the  nearly  acceleration-free  ('v  10-3  of 
Earth's  gravity)  environment  of  the  Spacelab.  In  such  surroundings  the  power  re- 
quirements relative  to  those  in  a terrestrial  laboratory  can  be  greatly  reduced  or 
the  size  range  of  controllable  drops  can  be  extended  (see  Section  2.3). 

The  utilization  of  radiation  pressure  as  a means  of  levitating  small 
spheres  has  been  successfully  demonstrated  in  the  laboratory  by  Ashkin  (1970). 

In  his  first  experiment,  an  argon  laser  with  an  output  of  a few  milliwatts  was 
used  to  accelerate  micron-sized  latex  spheres  which  were  freely  suspended  in  water. 
Ashkin  (1971)  later  demonstrated  the  stable  levitation  of  transparent  glass  spheres 
('v  20  pm  diameter)  in  air  and  in  vacuum  (^  1 Torr)  using  a laser  beam  of  250  mW. 
In  other  experiments  a 50  mW  laser  was  used  to  accelerate  'v  5 pm  diameter  water 
droplets  In  air. 
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During  his  laboratory  work,  Ashkin  (1970)  discovered  that  the  spheres  used 
in  his  experiments  were  "simultaneously  drawn  into  the  beam  axis  and  accelerated 
in  the  direction  of  light".  The  attraction  toward  the  center  of  the  beam  is  a re- 
sult of  the  gaussian  energy  distribution  across  the  beam  diameter  and  will  be  dis- 
cussed qualitatively  in  Section  2.1. A.  It  should  be  pointed  out  here  that  this 
effect  will  be  very  useful  in  maintaining  the  precise  position  of  droplets  since 
there  will  always  be  a force  tending  to  keep  the  droplet  in  the  beam  center  (sim- 
ilar to  a potential  well) . 

Utilization  of  optical  radiation  pressure  to  stably  suspend  small  parti- 
cles is  presently  being  considered  by  a group  at  the  Ruhr-Universitat  Bochum  in 
West  Germany  (Schwera,  1976).  Attempts  to  use  magnetic  forces  to  levitate  indivi- 
dual particles  for  experiments  on  Mie  scattering  were  complicated  by  the  oscilla- 
tion of  the  particles.  Experiments  are  now  being  designed  using  a laser  to  more 
stably  position  the  particles. 

The  capabilities  of  optical  radiation  pressure  as  a tool  for  studying  the 
microphysical  processes  involved  in  the  atmosphere  must  be  evaluated  in  terms  of 
the  experiments.  Typical  questions  to  be  answered  are: 

What  are  the  droplet  sizes  of  interest? 

To  what  velocities  must  the  droplets  be  accelerated? 

What  are  the  associated  laser  power  requirements? 

Are  these  power  requirements  compatible  with  the  constraints 

* 

of  the  laboratory  (in  the  case  of  Spacelab)? 

The  approach  of  this  study  will  be  to:  (1)  investigate  and  evaluate  the 
parameters  which  influence  optical  levitation  or  motion  control,  (2)  define  the 
limits  of  applicability  to  the  Atmospheric  Cloud  Physics  Laboratory  (ACPL) , (3) 
define  an  optical  motion  control  concept  for  ACPL,  and  (4)  outline  a program  for 
future  work  in  the  area  of  optical  motion  control. 

1 . 2 Parameters  Important  to  Optical  Levitation  or  Motion  Control 

The  efficiency  with  which  optical  radiation  pressure  can  be  applied  to 
the  levitation  and  positioning  of  particles  will  be  determined  by  the  character- 
istics of:  (1)  the  particles,  (2)  the  particle  environment,  and  (3)  the  laser. 

In  actuality,  these  factors  are  not  mutually  exclusive  but  for  the  purposes  of 
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the  following  discussion  their  relationship  to  the  present  objective  will  be 
treated  separately. 

I 

1.2.1  Particle  Characteristics 

The  particle  characteristics  which  obviously  affect  its  amenability  to 
levitation  include  its  mass  and  optical  properties  (i.e.,  reflection,  transmission 
and  absorption).  The  particle  mass  will  determine  the  applied  force  necessary  to 
suspend  or  accelerate  it  while  the  optical  characteristics  will  determine  the  ef- 
ficiency with  which  the  incident  optical  energy  is  converted  into  an  applied  force. 
It  should  be  noted  here  that  although  all  energy  absorbed  by  the  particle  contrib- 
utes directly  to  a momentum  transfer  to  the  particle,  this  absorbed  energy  also 
contributes  to  heating  of  the  particle.  Excessive  heating  of  the  particle  will, 
in  the  case  of  a water  droplet,  increase  its  evaporation  and  if  the  heating  results 
in  a large  temperature  gradient  in  the  particle  the  resulting  photophoretic  force 
may  become  comparable  to  the  applied  radiation  force.  The  importance  of  these 
considerations,  as  they  apply  to  optical  levitation  of  water  droplets,  will  be 
discussed  in  Section  2.1.2. 

1.2.2  Particle  Environment 

The  particle  environment  will  be  defined  to  consist  of:  (1)  the  medium  in 
which  the  particle  is  being  levitated,  and  (2)  the  ambient  acceleration  field 
(e.g.,  gravity)  which  is  present  in  the  laboratory.  The  temperature,  pressure  and 
viscosity  are  characteristics  of  the  medium  which  will  determine  the  thermophoretic 
forces  acting  on  the  particle  (Section  2.1.2)  as  well  as  the  maximum  velocity  a 
particle  can  attain  for  a given  applied  force  (Section  2.3.3).  The  opposing  accel- 
eration field  will  determine  the  laser  power  required  to  perform  the  actual  levi- 
tation. The  results  in  Section  2.3  illustrate  the  requirements  for  conditions  in 
a terrestrial  laboratory  as  well  as  for  a range  of  conditions  expected  to  exist  in 
the  Spacelab  on  Shuttle.  Figure  1 illustrates  the  sources  and  magnitudes  of  the 
perturbing  acceleration  fields  expected  for  Spacelab. 

1.2.3  Laser  Characteristics 

The  laser  characteristics,  other  than  power,  which  must  be  considered  in- 
clude the  efficiency  of  the  laser  and  the  optical  mode  of  the  beam.  The  laser 
efficiency  is  normally  not  the  controlling  factor  for  experiments  conducted  in  a 
terrestrial  laboratory  since  laser  input  power  and  physical  size  are  usually  not 
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FIGURE  1 . SHUTTLE/SPACELAB  ACCELERATION  LEVELS 


a problem.  However,  for  experiments  conducted  in  the  Spacelab,  power  and  size 
requirements  may  become  important.  Since  the  actua]  output  wavelength  has  little 
effect  on  the  required  levitation  power  (Section  2.3.3)  an  important  consideration 
in  choosing  a laser  (e.g..  Argon  vs  Helium-Neon)  should  be  its  efficiency. 

Selection  of  the  beam  mode  (TEM01  or  TEMqq)  is  dictated  largely  by  the 
requirements  of  the  experiment.  Askin  (1974)  has  successfully  levitated  spheres 
using  both  the  TEMqj  mode  (do-nut  shaped  beam)  and  the  TEMqo  mode  (gaussian  beam). 
His  data  shows  that  with  the  TEMqo  mode,  spherical  particles  can  be  levitated  much 
more  stably.  Since  precise  positioning  is  one  of  the  objectives  of  optical  motion 
control,  the  obvious  choice  for  consideration  is  the  TEMqq  mode  laser. 

1.2.4  Analysis  Method 

The  analyses  which  follow  evaluate  the  importance  of  the  above  parameters 
relevant  to  the  objective  of  levitating  (positioning)  and  controlling  the  motion 
of  water  droplets  in  air.  Consequently  the  detailed  examples  which  are  presented 
reflect  this  objective,  and  will  provide  guidelines  for  assessing  the  applicabil- 
ity of  optical  motion  control  to  problems  in  atmospheric  microphysics. 


2.0 


EVALUATION  OF  CRITICAL  FACTORS  FOR  OPTICAL  MOTION  CONTROL 


2 . 1 Forces  Relevant  to  Optical  Motion  Control 

2.1.1  Radiation  Pressure  Forces 

The  force  on  a body  due  to  radiation  pressure  (van  de  Hulst,  1957)  is 
given  by: 


F 


— Q 
c pr 


(1) 


where: 


and 


W = power  incident  on  the  body  (watts) 
c = speed  of  light  (cm  sec-*) 


pr 

The  value  of  Q 


pr 


efficiency  factor  for  radiation  pressure  (dimensionless) 
can  be  calculated  from  the  efficiency  factors  for  absorption 


Q , and  for  scattering  Q as: 

HDS  SC3 


Q = Q , + Q (1  - cos  0)  (2) 

xpr  xabs  xsca 

where  0 is  the  angle  between  the  scattered  light  and  the  transmitted  light  as 
illustrated  in  Figure  2.  Equation  (2)  indicates  that  all  absorbed  radiation  con- 
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FIGURE  2.  SCATTERING  ANGLE  DEFINITION 


tributes  directly  to  momentum  transfer,  while  only  the  backward  components  of 
scattered  radiation  (cos  0 = -1)  contribute.  However,  the  calculations  which 
follow  pertain  primarily  to  water  droplets  in  air  so  that  the  radiation  pressure 
due  to  the  Qa^s  term  is  negligible  in  comparison  to  the  other  terms,  i.e.,  water 
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droplets  are  considered  essentially  transparent  at  visible  wavelengths  (van  de 
Hulst,  1957).  Thus,  Equation  (2)  can  be  written  as: 


Q = Q (1  - cos  0)  (3) 

xpr  xsca 

where  Qgca  = Qext  (extinction  efficiency  factor)  for  a non-absorbing  sphere.  For 
the  case  of  water  droplets  with  index  of  refraction  m = 1.33  and  size  parameter 
x = 2-nr/A  >>  1,  van  de  Hulst  (1957)  gives  a value  of  Qext  = 2 (r  is  drop  radius 
and  A is  light  wavelength  in  same  units).  Equation  (3)  now  becomes: 

Qpr  = 2 (1  - cos  9)  (4) 

For  the  above  value  of  m,  van  de  Hulst  (1957)  gives  a value  of: 


(cos  6)  = 0.87 

which  when  substituted  in  Equation  (4)  gives  Qp^_  = 0.26  for  spherical  water  drop- 
lets. 

The  power  required  to  levitate,  i.e.,  exactly  counterbalance  the  gravita- 
tional force  on  a water  droplet  is  given  by  setting  F = mg  in  equation  (1)  which 
then  can  be  written  as: 


W 


_ mg  c 

Qpr 

4/3  it  p r3  gc 

V 


(5) 


which  gives  the  levitation  power  required  for  a particle  of  radius  r and  density  p 
in  the  earth's  gravitational  field.  For  a 10  pm  (radius)  spherical  water  droplet, 
a power  of  0.047  watts  incident  on  the  particle  surface  is  required  for  levitation. 
The  required  power  for  g = 980  cm/sec2  (at  earth's  surface)  is  plotted  as  a func- 
tion of  r in  Figure  3. 

Equation  (5)  can  be  converted  to  a power  density  (watts/cm2)  requirement 
by  merely  dividing  by  the  particle  area  to  obtain: 


1 


I = 


w 


(6) 


irr‘ 


= 4 p r gc 

3 Q 

Pr 

Thus,  the  beam  intensity  required  for  levitation  is  directly  proportional  to  the 
radius  of  the  droplet.  For  the  example  evaluated  above  (r  = 10-3  cm),  this  trans- 
lates into  a power  density  of  approximately  1.5  x 104  watts/cm2,  or  about  105 
times  the  intensity  of  solar  energy  on  the  earth  ('v  0.14  watts/cm2). 

It  should  be  noted  that  equations  (5)  and  (6)  are  general  and  can  be  ap- 
plied to  any  particle  for  which  the  proper  value  for  Q is  known  or  can  be  com- 
puted . 

2.1.2  Phoretic  Forces 

A temperature  difference  over  the  surface  of  a body  surrounded  by  a gas 
will  result  in  a net  force  on  the  particle.  If  the  uneven  heating  of  the  particle 
is  due  to  a temperature  gradient  within  the  gas  the  force  is  called  thermophoretic . 
If  a one-sided  illumination  of  the  particle  and  resulting  absorption  causes  an 
uneven  temperature  distribution  within  the  particle  the  resulting  force  is  called 
photophoretic. 

2. 1.2.1  Thermophoretic 

Consider  a sphere  residing  in  a gas  in  which  there  is  a temperature  gra- 
dient VT^  (e.g.,  in  a thermal  diffusion  chamber).  If  the  particle  radius  (r)  is 
much  greater  than  the  mean  free  paths  of  the  gas  molecules,  then  the  temperature 
gradient  within  the  spherical  body  is  given  by  Fuchs  (1964)  as: 


K 


3 k 

VT  = £ 


2 k +k(7Tg> 
g P 


(7) 


where:  k = thermal  conductivity  of  the  gas  (erg  cm-1  sec  1 °K-1) 

O 

kp  = thermal  conductivity  of  the  gas  (erg  cm-1  sec-1  °K-1 ) 
The  thermophoretic  force  on  the  particle  is  then  given  by  Fuchs  as: 
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7 I K 2 

F ^ £ — IL.  T—  (vtb)  (8) 

th  2k  + k P T ^ g;  K J 

8 P 8 8 

I 

where:  ri  = viscosity  of  the  gas  (gm  cm-1  sec-*) 

r = radius  of  spherical  particle  (cm) 

p = density  of  the  gas  (gm  cm-*) 

6 

T = temperature  of  the  gas  (°K) 

8 

The  calculated  value  for  F . below  illustrates  the  effect  of  the  thermo- 

th 

phoretic  force  on  a water  droplet  suspended  in  air.  The  following  values  were 
assumed  for  the  parameters  in  Equation  (8): 


k (air  at  21  °C)  = 
8 

kp  (H20  at  21  °C)  = 

n (air  at  21  °C)  = 

p (air  at  21  °C)  = 

8 

T = 294  °K 
8 


2.56  x 103  ergs  cm  1-sec-1-°K-1 
5.95  x 104  ergs  cm-1-sec-1-°K*'1 
1.838  x 10~4  gm  cm-1  sec-1 
1.201  x 10“3  gm  cm-3 


The  expression  for  F ^ then  becomes: 

F . = 1.07  x 10-7  r VT  (9) 

th  g 

Dividing  equation  (9)  by  the  mass  of  the  water  droplet  yields  an  expression  for 
the  thermophoretic  acceleration  given  by: 

1.07  x 10~7  r VT 

a = — (ic 

th  4/3  ir  p r3 


2.55  x 10“8 

r* 


Evaluating  the  above  expression  for  a 10  micron  water  droplet  radius  and  a temper- 
ature gradient  of  1 °K/cm  yields  a thermophoretic  acceleration  of 


11. 


= 2.55  x 10-?  cm  sec-2 
= 2.6  x 10-5  g / 


where  g is  the  earth's  gravitational  acceleration  of  the  980  cm/sec  . Equation 
(10)  indicates  that  the  temperature  gradients  in  the  medium  will  have  the  greatest 
effect  on  the  smallest  particles. 


The  Stokes  velocity  associated  with  this  acceleration  is  given  by  the 


formula : 


2 p r^ 

Vth  9 n 


which  gives  a value  of 


v , = 3. 1 x 10" b cm  sec' 
th 


2. 2. 2. 2 Photophoretics 

In  the  case  of  aspherical  particle  illuminated  from  one  side,  Fuchs 
gives  the  resulting  photophoretic  force  as: 

-3tt  q2  r R 

ph  2 PM  (k  + k ) i<x  * 

8 P 

where:  I = density  of  luminous  energy  flux  (erg  cm-2  sec-1) 

R = gas  constant  (ergs  °K-1  mole-1) 

0 

P = gas  pressure  (gm  cm-1  sec-2) 

M = molecular  weight  of  the  gas  (gm  mole-1) 
a = absorption  coefficient  of  the  particle  (cm-1) 

1 **  optical  absorption  thickness  (cm-1) 


and  all  other  parameters  are  as  previously  defined.  A typical  value  of  F is 

Ph 

calculated  by  assuming  the  following  values  for  the  parameters  in  Equation  (12) 
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= 8.317  x 107  ergs  °C-1  mole-1 


P = 1.10  x 106  gm  cm-1  sec-2 
M = 28.96  gm  mole-1 
a = 2.9  x 10“  3 cm-1 


The  computed  value  for  F ^ is  then: 


F , = -2.11  x 10-14  r I l gm-cm  sec 
ph 


-2 


(13) 


Dividing  equation  (13)  by  the  mass  of  the  water  droplet  gives  the  accel- 
eration due  to  the  photophoretic  force  as: 


ph 


-5.1  x 10-15 


1 *. 


(14) 


-1.0  x 10-14  - 
r 


where  a value  of  2r  (worst  case)  has  been  assumed  for  the  optical  absorption  thick- 
ness (l). 


Assuming  that  the  incident  levitation  power  W as  calculated  previously  is  focussed 
to  a beam  width  equal  to  the  particle  diameter  we  can  write: 


I = 


W 


(15) 


7rr^ 


and  Equation  (14)  becomes: 


a . = -1.0  x 10-14 


W 


ph 


(16) 


irrw 


Substituting  from  W from  Equation  (5)  gives: 

W: 


which  shows  that  a ^ is  independent  of  the  particle  size.  The  corresponding 
Stokes  velocity  for  a 10  p radius  water  droplet  in  air  has  a value  of: 


v , = -1.8  x 10-3  cm  sec-1 
ph 


In  the  above  computations  we  have  assumed  conditions  relevant  to  a water  droplet 
suspended  in  air. 

The  negative  sign  associated  with  the  photophoretic  force  and  acceleration 
indicate  that  the  particle  will  be  propelled  in  a direction  opposite  to  that  of 
the  incident  radiation  (negative  photophoresis)  whereas  thermophoretic  and  radia- 
tion pressures  are  in  the  direction  of  the  beam  for  water  in  air.  In  the  case  of 
nearly  transparent  particles  (such  as  water  droplets)  this  phenomenon  results  from 
the  fact  that  the  rays  refracted  inside  the  particle  will  cause  a greater  heating 
on  the  side  away  from  the  incident  light.  Highly  absorbing  particles,  on  the 
other  hand,  will  be  heated  more  on  the  front  surface  and  be  moved  in  the  direction 
of  the  incident  light  (positive  photophoresis) . 


2.1.3  Relative  Importance  of  Forces 


In  order  to  evaluate  the  influence  of  the  phoretic  forces  on  the  objective 
of  levitating  a particle  using  radiation  pressure,  the  ratios  apj1/apr  an<^  ath^apr 
can  be  calculated.  To  facilitate  this  computation  the  factors  upon  which  each 
acceleration  depends  are  presented  in  Table  I.  The  groupings  are  determined  by 
whether  the  parameters  are  characteristics  of  the  particle,  the  environment  or  the 
laser  (Section  1.2).  It  should  be  noted  that  a ^ depends  on  a^r  by  virtue  of  the 
fact  that  the  particle  heating  depends  directly  on  W.  Thus  it  can  be  concluded 
that  for  this  example  the  phoretic  forces  can  be  neglected  compared  to  the  radia- 
tion pressure  force.  Furthermore,  since  a^  varies  as  r-2  and  a^  is  independent 
of  r they  will  not  be  important  for  larger  droplets  (cf.  example  bottom  Table  I). 

2.1.4  Transverse  Stabilizing  Force  of  the  Laser  Beam 


The  forces  which  have  been  discussed  in  the  previous  sections  are  impor- 
tant in  terms  of  determining  whether  or  not  a particle  can  be  levitated  or  accel- 
erated using  optical  radiation.  However,  given  that  a laser  will  suspend  a par- 
ticle, why  does  the  particle  remain  in  the  beam  and  not  slowly  drift  to  the  side 
and  subsequently  fall?  As  previously  mentioned,  Ashkin  (1970)  observed  that  par- 
ticles were  actually  attracted  toward  the  center  of  the  laser  beam  (gaussian  mode) 


TABLE  I.  PARAMETERS  AFFECTING  VARIOUS  ACCELERATION  FORCES  ON  A PARTICLE 


T' 


A droplet  being  suspended  by  a laser  with  a gausslan  energy  distribution 
will  experience  a variation  in  the  incident  energy  across  its  surface.  Consider 
the  sketch  of  Figure  A which  illustrates  the  optical  paths  of  two  parallel  rays, 

Ij  and  I2,  of  light  incident  on  the  droplet.  Both  rays  are  equally  distant  from 
the  droplet  center  but  Ij  originates  nearer  the  center  of  the  beam  and  is  therefore 
stronger  than  I2.  The  forces  resulting  on  the  droplet  are  determined  by  drawing 
the  resulting  vectors  at  each  droplet-air  interface.  For  example,  the  incident 
ray  Ij  and  the  reflected  ray  Ijd  produce  a force  Fjm,  while  the  incident  ray  and 
the  refracted  ray  Ijr  produce  a force  Fuj^.  Similarly,  the  forces  resulting  at  the 
exit  point  of  ray  Ij  as  well  as  the  corresponding  forces  for  ray  I2  can  be  de- 
termined. 


k 

| 


As  shown  in  the  sketch,  all  forces  contribute  to  "pushing"  the  particle 
in  the  direction  of  the  incident  light.  However,  forces  Fjp^  and  Fid0  are  directed 
toward  the  center  of  the  beam  while  forces  F2Di  and  F2D0  are  directed  away  from  the 
beam  center.  The  forces,  FjRi  and  Fir0,  and  F2R1  and  F2R0  cancel  to  first  order 
(Ashkin,  1970).  Since  the  intensity  of  ray  Ij  is  greater  than  the  intensity  of 
ray  I2,  all  forces  associated  with  Ij  will  be  greater  than  the  corresponding  forces 
associated  with  I2.  Thus  we  can  write: 


or 


F > F 
IDi  2Di 


F > F 
lDo  2Do 


(18) 


F + F 
IDi  lDo 


F + F 
2Di  2Do 


(19) 


As  was  shown  in  Figure  A,  all  rays  incident  on  the  left  hand  of  the  drop- 
let will  result  in  a force  directed  toward  the  beam  center  (-x  direction)  while 
rays  incident  on  the  right  half  will  result  in  a force  directed  out  of  the  beam 
(+x  direction).  The  relative  magnitudes  of  all  rays  can  be  determined  by  their 
origin  within  the  gaussian  energy  distribution  of  the  beam.  For  example,  Ij  orig- 
inates at  the  beam  center  (x=0)  and  will  have  a magnitude  of : 


16. 


< 

: — >. 

FORCES  ON  LEFT  SIDE 

FORCES  ON  RIGHT  SIDE 

OF  DROPLET 

OF  DROPLET 

F1 Ri  + F1Ro  “ 0 

F2Ri  + F2Ro  = 0 

F1  Di  + F1Do  = fx- 

F2Di  + F2Do  = fx+ 

F = I f oc  A, 
X-  x-  1 

Fx+  = Z fx+  “ A2 

where  F = force  pushing  droplet  into  beam  center 

F^+  = force  pushing  droplet  out  of  beam  center 


Note:  Fx  /Fx+  « A1/A2  = Restoring  Force/Repelling  Force 

F^  - Fx+  oc  A1  - A2=  Net  Restoring  Force 

FIGURE  4.  ILLUSTRATION  OF  FORCES  ACTING  ON  A SPHERICAL 
DROPLET  BY  A GAUSSIAN  LASER  BEAM 


I,  “ — 
1 fa 


* 0.4 


while  I2  originates  at  the  2a  point  (x=l)  and  will  have  a value  of: 


^ 0.15 


However,  the  total  restoring  force  on  the  droplet  (F  ) will  be  proportional  to 
the  area  Aj  under  the  gaussian  energy  distribution  curve.  Similarly,  the  total 
force  tending  to  push  the  droplet  out  of  the  beam  (F  ) will  be  proportional  to  A2. 

The  relative  magnitudes  of  these  forces  can  easily  be  computed  for  various 
drop  sizes.  In  Figure  5a,  a droplet  of  diameter  la  (relative  to  the  width  of  the 
laser  beam)  is  positioned  in  a gaussian  laser  beam.  As  the  droplet  is  moved  across 
the  diameter  of  the  laser  beam,  the  corresponding  values  of  F^  /F^+  (i.e.,  Aj/A2) 
are  computed  in  relative  units  and  plotted  in  Figure  5b.  The  ratio  has  a minimum 
value  of  1 when  the  droplet  is  at  the  beam  center,  and  increases  as  the  droplet 
moves  to  either  side  of  center,  indicating  that  F^  is  always  greater  than  F^+. 
Similarly,  the  net  restoring  force  (-x  direction)  can  be  determined  qualitatively 
from  a plot  of  (F^_  - F^+)  as  shown  in  Figure  5c.  The  maximum  restoring  force 
occurs  when  the  droplet  is  centered  at  the  la  position  of  the  beam  corresponding 
to  the  point  at  which  the  intensity  gradient  of  the  gaussian  distribution  is  a 
maximum. 

Similar  results  are  plotted  in  Figure  6 for  a droplet  whose  diameter  is 
2a  (i.e.,  twice  as  large  as  the  above  example).  As  shown  in  that  figure,  which  is 
plotted  to  the  same  scale  as  Figure  5,  the  corresponding  values  of  (Fx_/Fx+)  anc* 

(F^  - F^+)  have  increased  greatly.  These  results  provide  qualitative  indication 

that  the  restoring  force  on  a particle  will  increase  as  the  droplet  size  and  beam 
size  become  comparable.  Figure  7 illustrates  this  effect  for  (F  - F ) , 

X ” Xt  D)3X 

i.e.,  a droplet  positioned  at  the  la  point  in  the  beam.  It  should  be  noted,  how- 
ever, that  the  restoring  force  will  become  less  effective  as  the  droplet  becomes 
so  large  that  it  effectively  presents  a flat  surface  to  the  beam. 
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6a.  DISTANCE  FROM  BEAM  CENTER  (UNITS  OF  a) 
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6b.  POSITION  OF  DROPLET  CENTER  RELATIVE  TO  BEAM 
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6c.  POSITION  OF  DROPLET  CENTER  RELATIVE  TO  BEAM 

FIGURE  6(a).  COMPARISON  OF  2a  WATER  DROPLET  WITH  GAUSSIAN  BEAM. 

(b) .  POSITIONAL  VARIATION  OF  RESTORING  FORCE/REPELLING 

FORCE  (Fx_/Fx+). 

(c) .  POSITIONAL  VARIATION  OF  NET  RESTORING  FORCE 

(px--Fx+>- 
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RELATIVE  UNITS 


12  3 4 

DROPLET  RADIUS  (UNITS  OF  a) 


FIGURE  7.  VARIATION  OF  MAXIMUM  RESTORING 
FORCE  AS  A FUNCTION  OF  DROPLET 
SIZE.  DROPLET  SIZE  IS  GIVEN  IN 
TERMS  OF  LASER  BEAM  WIDTH. 


The  foregoing  results  indicate  the  importance  of  minimizing  the  size  of 
the  laser  beam  as  it  impinges  upon  the  droplet.  In  addition  to  providing  the 
greatest  levitation  force,  it  will  also  maximize  the  "trapping"  effect  of  the 
gaussian  laser  beam.  This  last  effect  is  very  important  in  terms  of  stabilizing 
a body  and  will  enhance  the  use  of  optical  radiation  for  precise  positioning  of 
droplets. 

2.2  Heating  of  Optically  Levitated  Particles 

Consideration  must  be  given  to  the  heating  of  an  optically  levitated  par- 
ticle in  as  much  as  it  may  define  the  limits  within  which  optical  levitation  is 
feasible.  A fraction  of  the  incident  optical  energy  will  be  converted  into  ther- 
mal energy,  resulting  in  an  increase  of  the  droplet  temperature.  The  magnitude 
of  this  temperature  increase  will  depend  primarily  on  the  physical  characteristics 
of  the  droplet  and  the  medium  in  which  it  is  being  levitated . The  eventual  steady- 
state  temperature  can  be  computed  in  a relatively  straightforward  manner  by  equat- 
ing the  heat  input  (absorbed  optical  energy)  and  heat  loss  (radiation  and  conduc- 
tion to  the  medium)  for  the  levitation  power  required. 

The  heat  equivalent  of  the  incident  beam  can  be  expressed  by  the  equation: 
Q.n  = a k(l-q)  W (22) 

where:  = heat  input  to  the  droplet  (cal/sec) 

W = total  power  incident  on  droplet  (ergs/sec) 
q = fraction  of  incident  beam  which  is  reflected  (=  0.1) 

(1-q)  = fraction  of  incident  beam  transmitting  (=  0.9) 

k = conversion  factor  from  optical  energy  to  mechanical 
equivalent  of  heat  (=  0.239  x 10-7  cal/erg) 
a = fraction  of  transmitted  beam  which  is  absorbed  (2.9  x 10“  3 cm~*)Jl 
£.  = optical  absorption  thickness  (assume  = 2r,  diameter  of  drop 
as  a worst  case) 

Equation  (22)  can  thus  be  expressed  in  terms  of  the  levitation  power  and  the  drop- 
let size  as: 


Q,  = 1.25  x 10-10  Wr 
in 


(23) 


As  previously  stated  the  droplet  will,  given  sufficient  time,  attain  an 
equilibrium  temperature  at  which  the  heat  conduction  outward  from  the  center  of 
the  spherical  drop  can  be  expresses  as: 

Q = -k  A dT/dr  (24) 

out  B 

= -4 ir  r2  k dT/dr 
g 

which  will  be  true  for  any  drop  radius  (r).  In  the  above  expression  the  parame- 
ters are  defined  as  follows: 

• 

Q = rate  of  outward  heat  flow  (cal/sec) 
out 

r = radius  from  center  of  drop  (cm) 
kg  = thermal  conductivity  of  media  (erg  cm*'  sec-*  °K-*) 

T = absolute  temperature  (°K) 

A = spherical  surface  area  at  radius  r (cm2) 

Equation  (24)  can  be  integrated  to  obtain  the  outward  heat  flow  from  the  spherical 
surface  at  Temperature  Tg  to  the  conductive  medium  at  temperature  T with  the  re- 
sult: 


Q = 4n  r k (T  - T ) (25) 

c „ g s m7 

out 

Using  a value  of  6.1  x 10-*7  cal  cm-*  sec-* °K-*  for  the  thermal  conductivity  of  air 
we  obtain: 

Q = 7.7  x 10-4  r (T  - T ) (26) 

c s m 

out 

Similarly,  the  radiative  heat  loss  from  the  droplet  can  be  calculated  in 
a straightforward  manner  using  the  equation: 

% • ' 0 <Ts  - ■£>  (27> 

out 

where:  e = emissivity  of  the  water  droplet 

a = Stephen-Boltzman  constant  (1.36  x 10-*2  cal  sec-*  cm-2  °K-4 ) 


A 
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The  emissivity  is  given  by: 

e = 1 - e-a*  (28) 

^ a H for  a Si  « 1 

where:  a = absorption  coefficient  of  the  drop  (2.9  x 10-3  cm-1  i.) 

and  Jt  = optical  absorption  thickness  (assume  t « 2r,  diameter 

of  drop  as  a worst  case) 

Substitution  into  Equation  (6)  yields: 

Q = 7.9  x 10"15  r (T4  - T4)  (29) 

xr  s m 

out 

The  relative  importance  of  radiative  and  conductive  heat  loss  from  the 
droplet  is  readily  obtained  by  dividing  Equation  (25)  into  Equation  (29)  to  ob- 
tain: 


IQ-n  (t3  + T2  T + T T2  + T3)  (30) 

s s m s m m 

Assuming  a worst  case  condition  of  Tg  = 373  °K  (boiling  point  of  water)  and  T^  = 

294  °K  (room  temperature)  we  find: 

^ 1.5  x 10~3. 

Clearly,  the  radiative  heat  loss  is  insignificant  when  compared  to  the  conductive 
heat  loss  to  the  surrounding  medium  (assumed  to  be  air  at  1 atm).  However,  it 
must  be  pointed  out  that  if  the  levitation  is  attempted  at  successively  lower  pres- 
sures the  conductive  effects  will  decrease  until  eventually  the  radiative  term  be- 
comes dominant. 

As  a consequence  of  Equation  (30)  the  equilibrium  temperature  of  the  drop- 
let can  be  calculated  by  considering  the  simplifed  heat  balance  equation: 
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(31) 


«in  " <>c  , 
out 


or 


1.25  x 10~10  Wr  = 7.7  x 10~4  r (T  - T ) 

s m 


(32) 


which  is  obtained  from  Equations  (23)  and  (26).  Solving  for  the  temperature  in- 
crease of  the  drop  we  have: 


AT  = (T  - T ) = 1.6  x 1 0-7  W (33) 

s m 

where  W (the  levitation  power)  is  given  by  Equation  (5)  as: 

w = 4/3  it  p r3  sc 


The  temperature  delta  then  becomes: 

AT  = 7.7  x 104  g r3 


(34) 


where  a value  of  p = 1 (water)  has  been  assumed. 

Thus,  for  a 10  ym  (radius)  spherical  water  droplet,  being  levitated  at 
standard  temperature  and  pressure  conditions,  in  a terrestial  laboratory,  a tem- 
perature increase  of  approximately  0.08  °C  will  be  observed.  A graph  of  required 
levitation  power  and  the  corresponding  temperature  rise  is  given  in  Figure  8 as  a 
function  of  particle  size  assuming  opposing  acceleration  fields  of  1 g,  10-3  g and 
10-6  g,  where  g = 980  cm/sec2. 

2.3  Practical  Limits  of  Acceleration  and  Velocity 
2.3.1  Typical  Acceleration  Levels 

As  discussed  in  Section  2.1.3,  the  acceleration  which  can  be  imparted  to 
a particle  by  optical  radiation  pressure  is  expressed  by  the  equation: 


a 


pr 


3 

4ctt 


ir 

pr3 


W. 


(35) 
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POWER  REQUIRED  FOR  LEVITATION  (WATTS) 
INCREASE  IN  DROP  TEMPERATURE  ABOVE  AMBIENT 


Substituting  the  values  of  Q - 0.26  and  p = 1 gm/cm3,  which  are  relevant  to  a 
water  droplet,  Equation  (35)  becomes: 

a = 2.07  x 10-12  W/r3.  (36) 

pr 


The  variation  of  a with  changing  values  of  r and  W can  now  be  easily  calculated 
by  alternately  holding  r and  W at  constant  values.  Figure  9 illustrates  the  vari- 
ation of  radiation  pressure  acceleration  (a  ) as  a function  of  droplet  radius  for 
several  values  of  laser  power  input.  Similarly,  Figure  10  indicates  a as  a func- 
tion of  W for  several  values  of  r. 

In  order  to  obtain  the  net  acceleration  of  the  particle,  the  opposing  ac- 
celeration (e.g.,  earth's  gravity)  must  be  subtracted  from  the  values  plotted  in 
the  figures.  For  example,  working  in  a terrestrial  laboratory  with  a 100  milli- 
watt laser.  Figure  9 indicates  that  an  acceleration  of  1 g (980  cm/sec2)  can  be 
given  to  a 14  pm  radius  droplet.  However,  this  is  just  equal  to  the  opposing  ac- 
celeration of  earth's  gravity  so  that  the  net  acceleration  of  the  droplet  is  zero. 
If  the  same  laser  were  used  in  the  Spacelab  where  a typical  acceleration  level 
might  be  *v<  10“  3 g (Figure  1)  the  net  acceleration  of  the  particle  would  be  almost 
1 g.  Similarly,  Figures  9 and  10  can  be  used  to  determine  the  net  acceleration  on 
a water  droplet  for  various  values  of  droplet  size,  laser  input  power  and  opposing 
acceleration  field. 


2.3.2  Wavelength  Dependence  of  Levitation  Power 


The  particle  acceleration  achieved  is  primarily  a function  of  the  Mie  size 
parameter,  x = 2ir  r/A,  which  relates  the  particle  radius  (r)  and  the  wavelength 
(A).  Maximum  acceleration  is  achieved  with  x = 1 (Adler,  et  al.,  1972)  with 
smaller  local  maximum  at  multiples  of  this  value.  In  the  present  case  of  optical 

o 

levitation  (A  = 6328  A)  being  considered  herein,  x >>  1 for  particles  greater  than 
one  micron.  As  a result,  the  value  of  Q = 0.26  will  remain  virtually  unaffected 
as  A is  varied  within  the  bounds  of  the  visible  spectrum,  and  the  levitation  power 
required  will  remain  constant. 


For  completeness,  the  acceleration  magnitudes  achievable  in  a vacuum  from 
a broad  spectrum  of  electromagnetic  radiation  sources  is  presented  in  Table  II 
(Adler,  et  al. , 1972).  In  all  cases  the  power  density  was  assumed  to  be  1 watt/cm2 
and  the  particle  density  1 gm/cm3.  As  a comparison,  acoustical  motion  control  as 
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FIGURE  9.  ACCELERATION  DUE  TO  RADIATION  PRESSURE  AS  A 
FUNCTION  OF  WATER  DROPLET  SIZE  FOR  SEVERAL 
LASER  POWER  LEVELS. 
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TABLE  II.  ACCELERATION  MAGNITUDES  AS  A FUNCTION  OF  WAVELENGTH  OF  INCIDENT  ENERGY 
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being  developed  by  Jet  Propulsion  Laboratory  has  effective  wavelengths  of  greater 
tha  3 cm  and  therefore  is  restricted  to  particle  radii  greater  than  a centimeter. 
From  this  table,  infrared  and  visible  electromagnetic  radiation  are  the  most  suit- 
able for  motion  control  of  the  micrometer  radii  of  interest  in  atmospheric  precip- 
itation processes  involving  the  growth  transition  zone  between  pure  diffusional 
growth  and  the  dynamic  collision/coalescence  growth. 

Thus,  along  with  the  ability  to  operate  at  a distance,  manipulate  selected 
particles  in  a group  of  particles  and  minimal  effects  on  the  particle  physical 
characteristics,  the  visible  and  near  visible  radiation  exerts  the  optimum  forces 
on  the  particles  of  interest  as  well  as  being  available  at  the  power  densities 
(lasers)  required  to  be  implemented. 


2.3.3  Typical  Velocity  Levels 

A water  drop  being  accelerated  in  air  by  a constant  applied  force  will 
eventually  attain  an  equilibrium  velocity  (v)  determined  by  the  resistive  drag 
force  which  can  be  expressed  as  follows  (Mason,  1971): 
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(37) 


where  is  the  drag  coefficient  and  is  a function  of  the  Reynolds  number  Re> 
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(38) 


defined  by  the  equation: 


C,  R Fa 

d e d 
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(39) 


The  drag  force  can  then  be  expressed  as 


/C,  R i 

Fd  = -§^|  6 * n r v 


(40) 


and  the  force  equation  for  the  drop  becomes: 


li 
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(41) 


The  left  side  of  Equation  (41)  is  just  the  net  force  acting  on  the  droplet.  The 
first  term  on  the  right  is  the  accelerating  force  acting  on  the  droplet  (neglect- 
ing the  buoyancy  of  air),  while  the  second  term  in  the  previously  defined  drag 
force,  at  equilibrium  (i.e.,  dv/dt  = 0).  Equation  (41)  can  be  solved  for  the 
velocity  to  obtain: 
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(42) 


Substituting  for  the  particle  mass  we  have: 


v = 


24 

2 pr2 

C , R 

9 n 

d 

e/ 

(43) 


In  the  case  of  Stokes  Law  (Rg  - 0.1),  R^  = 24  and  Equation  (39)  becomes: 


v = a P 


2 r2a 


9 K n 


(44) 


As  R^  becomes  larger,  Stokes  Law  increasingly  overestimates  the  actual  velocity  so 
that  empirical  relations  between  R^  and  R^  are  usually  used.  The  data  of  Beard 
and  Pruppacher  (1969)  showed  that  water  drops  can  be  well  represented  by  the  fol- 
lowing formulae  given  by  Mason  (1971) : 


C R 

■ v;-  -■  = 1 + 0.10  R0-955,  0.01  < r < 2,  (45a) 

24  e e 

= 1 + 0.11  R°-81,  2 - Rg  - 21,  (45b) 

= 1 + 0.189  R°'632,  21  - R - 200 
e e 


In  practice  R2  is  usually  plotted  as  a function  of  R^  as  shown  in  Fig- 
ure 11.  The  quantity  of  R2  can  be  expressed  in  terms  of  the  drop  radius,  the 
air  (ambient)  characteristics  and  the  acceleration  level  as  (Beard  and  Pruppacher, 
1969): 


C R2  = r3  p p —2 
d e 3 g n" 


(46) 
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where  a is  the  net  acceleration  of  the  particle.  Thus  the  procedure  for  calcu- 
lating the  equilibrium  velocity  of  a water  droplet  is  as  follows: 

1.  Specify  the  drop  characteristics  (r,p),  the  ambient 

characteristics  (p  ,n),  and  the  net  acceleration  (a). 

8 

2.  Calculate  C,  R2  from  Equation  (46) . 

a e 

3.  Determine  the  corresponding  value  of  Rg  from  Figure  11. 

4.  If  Re  < 0.1  calculate  v from  Equation  (44). 

5.  If  R - 0.1  calculate  v from  Equation  (45). 

e 

Figure  12  illustrates  the  equilibrium  or  terminal  velocity  attained  by  a water 
droplet  in  air  (T  = 294  °K,  P = 1 atm)  for  various  multiples  of  earth's  gravity 
(g). 

In  order  to  utilize  laser  energy  to  control  the  motion  and  position  of  a 
droplet,  the  acceleration  imparted  to  it  must  be  sufficient  to  exceed  the  opposing 
acceleration  field  (e.g.,  earth's  gravity  in  a terrestrial  laboratory).  In  Figures 
13,  14,  and  15  the  equilibrium  velocity  has  been  plotted  as  a function  of  droplet 
radius  for  various  values  of  laser  input  power.  The  three  figures  illustrate  the 
limitations  imposed  by  three  different  opposing  acceleration  fields,  i.e.,  1 g, 

10-3  g,  and  10-6  g,  respectively.  The  second  and  third  acceleration  levels  are 

typical  of  what  may  be  called  the  random  and  known  values,  respectively  (refer  to 
Figure  1) . The  surrounding  medium  is  again  assumed  to  be  air  at  T = 294  °K  and 
P = 1 atm. 

In  each  figure  the  equilibrium  velocities  which  would  pertain  if  there 
were  no  opposing  acceleration  are  indicated  by  the  dashed  lines.  Also  shown  on 
these  figures  is  the  terminal  velocity  a droplet  would  have  in  the  earth's  atmo- 
sphere for  the  conditions  of  1 g,  T = 294  °K  and  P = 1 atm.  Thus,  if  it  is  de- 
sired to  simulate  droplet  collisions  resembling  those  which  might  be  naturally 
occuring  in  the  atmosphere,  the  "1  g"  dot-dash  line  shown  in  the  figures  can  be 
used  as  a guide  in  determining  the  laser  power  required  for  a given  drop  size.  If 
a particular  experiment  requires  precise  positioning  of  a droplet,  the  droplet  ve- 
locity must  be  relatively  low,  e.g.,  < 0.1  cm /sec.  The  optimum  laser  power  to 
achieve  this  goal  can  be  obtained  by  referring  to  the  proper  figure. 
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DROPLET  VELOCITY  (CM/SEC) 


FIGURE  14.  WATER  DROPLET  VELOCITY  AS  A FUNCTION  OF  DROPLET 
SIZE  (OPPOSING  ACCELERATION  = 1 0“3g ) . 
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FIGURE  15.  WATER  DROPLET  VELOCITY  AS  A FUNCTION  OF  DROPLET 
SIZE  (OPPOSING  ACCELERATION  = 1 0“6g ) . 


Similar  results  are  presented  in  Figures  16,  17,  and  18.  However,  in 
this  case  the  independent  variable  is  the  laser  input  power  and  the  equilibrium 
velocity  is  presented  for  a range  of  droplet  sizes.  These  figures  can  be  used  by 
themselves,  or  in  combination  with  the  three  previously  discussed  figures  to  de- 
fine the  laser  input  power  appropriate  for  a given  experiment.  The  obvious  point 
illustrated  by  these  figures  is  however  the  increased  versatility  and  applicabil- 
ity of  optical  motion  control  as  one  moves  from  a terrestrial  laboratory  to  Space- 
lab. 
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LASER  INPUT  POWER  (WATTS) 


WATER  DROPLET  VELOCITY  AS  A FUNCTION  OF  LASER 
POWER  (OPPOSING  ACCELERATION  = 10"3g). 


FIGURE  17 
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APPLICABILITY  OF  OPTICAL  MOTION  CONTROL  TO  CLOUD  PHYSICS  STUDIES 


3 . 1 Atmospheric  Cloud  Physics  Laboratory  on  Shuttle/Spacelab. 

The  previous  sections  have  shown  that  the  utilization  of  optical  motion 
control  (OMC)  is  primarily  determined  by  the  available  laser  power.  The  Air  Force 
(among  others)  is  currently  developing  a space  qualified  laser  (Program  405B)  for 
use  in  satellite-to-satellite  communications  which  will  be  compatible  for  use  on 
Shuttle/Spacelab.  This  laser  will  operate  at  a wavelength  of  0.532  pm  with  a power 
output  of  between  100  to  200  raW  and  will  be  available  around  1979.  We  can,  there- 
fore, assume  that  the  first  application  of  OMC  to  the  ACPL  will  have  a useable 
laser  power  of  100  mW.  In  addition,  we  can  assume  that  the  maximum  residual  or 
ambient  acceleration  in  the  Spacelab  is  10-3  g (cf.  Figure  1).  Within  these  bounds 
we  can  define  the  applicability  of  OMC  for  experiments  relating  to:  (1)  the  growth 
of  water  droplets  or  ice  crystals,  and  (2)  the  collision  or  coalesence  of  two  or 
more  drops  or  crystals. 

The  role  of  OMC  in  the  first  of  these  applications  will  depend  on  the  meth- 
od by  which  the  droplets  or  crystals  are  generated.  For  example,  they  may  be 
generated  "in  position"  within  the  supersaturated  environment  of  a static  diffusion 
chamber.  The  laser  beam  would  then  be  used  as  a probe  to  select  and  hold  an  indi- 
vidual droplet  in  position  so  that  its  growth  history  could  be  observed.  An  alter- 
nate method  would  be  to  generate  droplets  singly  (e.g.,  using  a wire  probe  retrac- 
tor) near  the  edge  of  the  chamber  and  then  use  the  laser  energy  to  move  the  droplet 
into  a location  where  it  could  be  fixed  in  a viewing  position  by  a second  laser 
beam.  The  droplet  size  range  to  which  these  applications  can  be  applied  can  be 
easily  visualized  by  referring  to  Figure  19  (this  is  a modified  version  of  Figure 
14,  Section  2.3.3). 

Droplet  sizes  in  the  range  of  5-100  pm  radius  are  of  particular  interest 
for  cloud  physics  investigations  on  the  microphysical  level.  The  cross-hatched 
area  in  Figure  19  indicates  the  maximum  velocities  which  can  be  expected  for  an 
applied  laser  power  of  100  mW.  Droplets  with  radii  less  than  % 13  pm  can  be  accel- 
erated to  velocities  which  are  greater  than  their  respective  terminal  earth's  atmo- 
sphere (1  g)  velocities.  Dynamic  processes  such  as  scavenging,  and  collision/co- 
alescence can  therefore  be  simulated  and  observed  in  zero  g as  they  might  naturally 
occur  in  the  1 g earth  environment  but  with  better  control  and  longer  observation 
times.  Figure  19  also  illustrates  that  the  upper  size  limit  for  position  control 


43 


DROPLET  VELOCITY  (CM/SEC) 


MOTION 

CONTROL 


FIGURE  19.  WATER  DROPLET  VELOCITV'AS  A FUNCTION  OF  DROPLET 
SIZE.  CROSS-HATCHED  AREA  INDICATES  REGION  OF 
APPLICABILITY  FOR  SHUTTLE/SPACELAB  EXPERIMENTS 
(ASSUMING  100  mW  LASER). 
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of  water  droplets  will  occur  at  'v  130  urn  (the  size  for  which  the  100  mW  curve  be- 
comes essentially  vertical).  This  size  limit  could  also  be  obtained  from  Figure  9 
by  noting  the  droplet  size  at  which  its  radiation  pressure  acceleration  is  equal 
to  the  opposing  ambient  acceleration  of  10~3  g. 

Thus  OMC  can  provide  a significant  contribution  to  atmospheric  microphysics 
experimentation  in  ACPL  through  positioning  and  variable  accelerated  motion  con- 
trol. 

3 . 2 Terrestrial  Laboratory 

In  principle,  the  utilization  of  OMC  to  perform  cloud  physics  experiments 
in  a terrestrial  laboratory  is  not  limited  by  the  available  laser  power.  However, 
as  was  shown  in  Figure  8,  the  droplet  temperature  increases  rapidly  as  the  laser 
input  power  is  increased.  Thus,  the  applicability  of  OMC  will  be  determined  by 
the  acceptable  temperature  rise  of  the  droplet.  For  example,  if  the  maximum  tol- 
erable droplet  temperature  increase  were  10  °C  then  the  maximum  droplet  size  cutoff 
would  occur  at  'v  50  pm  radius  (laser  power  imput  of  'v  700  mW) , which  is  well  below 
the  previously  discussed  limit  of  'v  130  pm  for  the  ACPL  on  Shuttle.  Note  that  from 
Figure  8 the  power  required  to  levitate  a 100  pm  water  droplet  would  also  raise  its 
temperature  to  the  boiling  point  (1  g conditions). 

On  the  basis  of  the  above  considerations  terrestrial  applications  of  OMC  to 
water  droplets  will  be  confined  to  radii  less  than  'v  50  pm.  However,  earth  based 
laboratory  experiments  with  OMC  will  be  necessary  to  provide  proper  development  of 
motion  and  position  control  techniques  which  may  subsequently  be  implemented  on  the 
Shuttle/Spacelab  ACPL. 
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4.0  OPTICAL  MOTION  CONTROL  CONCEPT  FOR  ACPL 

A summary  of  potential  ACPL  experiments  which  would  benefit  from  the  imple- 
mentation of  optical  motion  control  capability  is  presented  in  Table  111.  A poten- 
tial experiment  illustrating  a specific  application  of  OMC  is  outlined  below. 

Assume  that  the  experiment  objective  is  to  monitor  the  growth  history  of 
an  individual  water  dropler  (or  ice  crystal)  for  an  extended  period  of  time.  A 
typical  protocol  would  be  as  follows: 

1.  A 25  pm  radius  water  droplet  is  generated  at  the  outer  edge 
of  a chamber  (e.g.,  static  diffusion  chamber  ) which  has  a 
radius  of  15  cm. 

2.  A 100  mW  laser  is  then  used  to  move  the  droplet  to  the  center 
of  the  chamber,  a distance  of  15  cm,  in  approximately  15  sec 
(from  Figure  19). 

3.  The  droplet  is  then  held  in  this  position  by  an  additional 
laser  beam  of  equal  power  so  that  its  growth  history  can  be 
monitored.  (Note:  the  diffraction  pattern  from  the  droplet 
can  be  used  for  a precise  measurement  of  its  size.) 

4.  The  droplet  will  continue  to  grow  until  its  radius  is  approx- 
imately 130  pm.  As  it  approaches  and  exceeds  this  size,  am- 
bient acceleration  forces  will  eventually  dominate  the  applied 
radiation  pressure  forces  and  the  droplet  will  be  lost. 

A sketch  illustrating  the  experiment  configuration  is  given  in  Figure  20. 

Similar  scenarios  can  be  envisioned  for  droplet  collision/coalescence 
experiments  where  impact  parameters  can  be  predetermined  to  better  than  one  microm- 
eter. The  droplets  can  be  accelerated  toward  each  other  or  in  the  same  direction 
with  velocities  representative  of  lg,  0.1  g,  etc.  acceleration  fields.  The  repeat- 
able  and  precise  control  will  permit  a better  assessment  of  the  effects  of  charge, 
surfactants,  etc.  on  droplet/droplet  and  droplet/ice  interactions. 
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TABLE  III.  ACPL  EXPERIMENT  MOTION  CONTROL  REQUIREMENTS 


Class  20  Unventilated  Droplet  Position  control  for  extended  observation 

Diffusion  Coefficient  periods 


r • 


MICROSCOPE  FOR 
DROPLET  OBSERVATION 


L - Focussing  lens 


• Stability  in  x direction  provided  by  opposing  laser  beams 

• Stability  in  y and  z direction  provided  by  transverse  restoring 
force  of  gaussian  beam 


FIGURE  20.  ILLUSTRATION  OF  THE  EXPERIMENT  CONFIGURATION  FOR 
GENERATING  AND  POSITIONING  A WATER  DROPLET  TO 
DETERMINE  ITS  GROWTH  HISTORY. 
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5.0 


PROGRAM  FOR  FUTURE  WORK 


This  study  indicates  that  continued  investigation  (primarily  in  the  lab- 
oratory) of  the  potential  application  of  OMC  to  atmospheric  cloud  physics  problems 
will  provide  valuable  insight  into  several  microphysical  phenomena.  For  example, 
experimental  determination  of  the  coalescence  efficiencies  of  cloud  droplets  in  the 
growth  transition  region  between  8 and  15  ym  radius  has  been  severely  hampered  due 
to  the  lack  of  sufficiently  accurate  control  of  these  small  droplets.  As  discussed 
in  the  previous  section  these  sizes  are  well  within  the  range  for  which  OMC  is  ap- 
plicable to  a terrestrial  laboratory.  In  addition.  Dr.  Ashkin  of  the  Bell  Tele- 
phone Laboratories  has  demonstrated  the  precise  positioning  capability  (positioning 
accuracies  of  a few  micrometers)  of  OMC  in  his  laboratory.  The  use  of  optical  lev- 
itation will  permit  the  establishment  of  accurate  initial  relative  positions  between 
two  drops  which,  upon  release,  will  undergo  a collision/coalescence  process.  The 
benefits  of  this  technique  are  further  enhanced  when  it  is  used  in  conjunction  with 
a drop  generator  (such  as  a wire  probe  retractor)  which  can  repeatedly  generate 
individual  drops  with  a given  size  and  trajectory  (i.e.,  fall  velocity).  It  will 
then  be  possible  to  precisely  control  the  impact  parameter  permitting  the  investi- 
gation of  various  droplet  interactions  as  the  impact  parameter  is  varied.  Thus, 
the  influence  of  droplet  charge,  electric  fields,  surfactants  and  turbulence  on  the 
collision/coalescence  process  in  the  size  range  between  8 and  15  ym  radius  can  be 
studied.  The  investigations  will  also  provide  a means  of  evaluating  the  utility  of 
optical  motion  control  for  experiments  in  the  ACPL.  In  addition,  this  levitation 
technique  may  be  extended  to  the  study  of  ice  crystal  electrification  and  coagula- 
tion (sticking  coefficient  determination). 

The  basic  equipment  to  implement  these  concepts  currently  exists  in  the 
General  Electric  Space  Sciences  Laboratory.  In  particular  the  following  items  are 
available: 

1.  50  mW  Helium-Neon  laser. 

2.  2.5  W Argon  laser. 


Aerosol  generators  and  atomizers  for  producing  a cloud  of  droplets. 

Wire  probe  retractor  and  vibrating  orfice  generators  for 
generating  individual  drops. 

Microscopes  for  observing  droplet  interactions. 
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A program  for  laboratory  evaluation  of  the  OMC  concept  as  applied  to  atmo- 
spheric cloud  physics  is  given  below. 

Experimental  Objectives 

1.  Generate  a cloud  of  droplets  and  demonstrate  the  capability  of 
capturing  and  levitating  individual  droplets  composed  of  a non- 
volatile substance  such  as  silicon  oil. 

2.  Repeat  objective  1 using  water  droplets.  Thit  objective  will 
require  the  design  of  a container  to  provide  a saturated 
environment  for  the  droplet  to  prevent  evaporation. 

3.  Implement  the  wire  probe  reactor  to  operate  in  conjunction 
with  the  optical  levitation  technique  (I.e.,  generate  and 
capture  individual  drops  with  repeatability) . 

4.  Utilize  the  wire  probe  retractor  to  generate  drops  which  fall 
past  another  drop  held  in  position  by  the  laser  beam.  The  ef- 
fect of  varying  drop  size  and  impact  parameter  can  be  studied. 

5.  Utilize  the  procedure  of  objective  4 to  investigate  the  in- 
fluence of  droplet  charge,  electric  fields,  surfactants  and 
turbulence  on  the  collision/coalescence  process. 

6.  Implement  a technique  recer  ,±y  developed  by  Ashkin  and 
Dziedzic  (personal  communication,  1977)  to  provide  feedback 
stabilization  of  optically  levitated  particles. 

It  should  be  noted  that  the  above  program  objectives  are  structured  in  a 
manner  which  permits  them  to  be  accomplished  incrementally,  i.e.,  each  objective 
contributes  to  those  which  follow.  The  exception  is  objective  6 which  is  rela- 
tively independent  of  the  others.  This  technique  is  most  important  when  working 
at  low  pressures  where  the  viscous  damping  forces  become  small  and  for  extended 
operation  in  zero  g where  the  laser  power  is  turned  up  periodically  for  reposi- 
tioning as  needed  and  turned  to  a low  level  standby  to  provide  position  sensing. 

In  conclusion,  useful  scientific  information  is  provided  by  each  of  the  first  five 
objectives  while  also  contributing  toward  evaluating  the  OMC  concept  for  use  in 
the  ACPL  on  Shuttle/Spacelab. 
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